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HIGH OBSERVED BRIGHTNESSES IN RADIO JETS 


M. SPADA^ M. SALVATP AND F. PACINT’^ 

ABSTRACT 

We study the variability properties at radio frequencies of the jets thought to be 
typical of Active Galactic Nuclei, i.e. bulk Lorentz factor P ~10, and incoherent syn¬ 
chrotron emission. We assume that the radiating electrons are accelerated at shocks 
within the jet, and that these shocks have a suitable (conical) geometry. In our frame¬ 
work, we can reproduce the variability pattern of Intra Day Variables (20% in less than 
1 day) as long as the observed brightness temperature is Tg < 3 10^^ K. The only strin¬ 
gent condition is an injection timescale of the perturbation shorter than the variability 
timescale; the geometric condition on the viewing direction is not especially severe, 
and agrees with the observed occurrence of the phenomenon. For higher values of Te 
coherent processes are probably necessary. 


Subject headings: BL Lacertae objects : individual (Mkn 421, S5 0716-1-714) — galaxies 
: jets — radio continuum : galaxies — radiation mechanisms : non thermal 


1. INTRODUCTION 

One of the defining features of Active Galactic Nuclei (AGN’s) is flux variability over most of 
the electromagnetic spectrum, with observed timescales which range from minutes up to several 
years ([Ulrich, Maraschi &: Urry 1997 ); the information content of the variations makes them a 
fundamental clue for the understanding of AGN’s. In a previous paper ( ^alvati, Spada &: Pacini 
|1998[ hereafter Paper I) we have considered the extremely fast variations (duration of about 20 
minutes) observed in the TeV range from Mkn 421 ( jGaidos et al. 1996 ). These flares imply an 
extremely compact emitting region and a very high photon-photon opacity, unless the source is 
moving with bulk Lorentz factor F ~ 100 (see also [Gelotti, Fabian &: Rees 199^ ). We have shown 
that extreme jet parameters are not needed, however, if the photon distribution is anisotropic in 
the comoving frame; such a distribution is a natural consequence of conical shocks, with opening 
angle « l/F. 

At the opposite end of the e.m. spectrum, similar problems are presented by fast radio vari¬ 
ations. Although radio emission occurs at much larger scales than TeV emission (10^® — 10^^ cm 
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versus 10^® — 10^'^ cm), variability up to a level ~ 20% has been observed in various sources with 
timescales shorter than a day, the so called Intra Day Variability (IDV) ([Wagner &: Witzel 1995 


hereafter WW). The observations imply brightness temperatures up to 10^® —10^® K, well above the 
Compton limit 10^^ K ( Kellermann Sz Pauliny-Toth 196^ ). Because of the relativistic motion of the 
source the intrinsic Tg is lower than the brightness temperature derived from the light curves 
(WW; [Rees 19^ ). Again, bulk Lorentz factors T « 100 would be required to cure the problem 
(see also Begelman, Rees Sz Sikora 1994[ ). VLBI observations indicate T ~ 10 only (jVermeulen fc 
[Cohen 1994 ; Ghisellini et al. 1993 ), so the intrinsic Tg remains well above 10^^ K, even after the 
beaming effects have been corrected for. 

Our purpose in this paper is to show that the model developed earlier to interpret the high 
energy variability may be relevant to the radio variability as well, with some appropriate modifi¬ 
cations. We shall show that in this way one can account for Te up to several times 10^^ K. Higher 
values of Tb wich have been claimed would however defy our interpretation. This is especially 
true for PKS0405-385 ( Kedziora-Chudczer et al. 1997| ), which varied on timescales of minutes with 
Tb ~ 10^^ K. Barring a relativistic jet with T ~ 10^, here one is compelled to accept intrinsically 
high brightnesses, due perhaps to coherent processes. 


2. DESCRIPTION OF THE MODEL 

Light travel time effects are important in the modelling of variability since they change the 
observed time with respect to the coordinate time by the factor 1 — f5cos9, with /3 the velocity 
of the source in units of c, and 9 the angle between the velocity vector and the line of sight. If 
the relevant source velocity is a physical velocity, one recovers the well known “shortening” factor 
~ 1/r^. In our scheme, the perturbation is modelled as a slab of extra electrons superposed on 
the steady jet, and flowing down the jet with the same bulk Lorentz factor T; the extra electrons 
radiate only after having passed through an oblique (conical) shock, where they are accelerated 
to relativistic energies in the comoving frame; that is, instead of radiating simultaneously, they 
radiate in succession according to their distance from the shock apex (see Fig. 1 for details). The 
relevant source velocity, then, is the phase velocity of the intersection point between the slab and 
the shock; calling a the incidence angle, we rewrite the “shortening” factor as 1 —/3cos(0 —a)/cosa, 
with P the physical velocity corresponding to T, and P/cosa the phase velocity. The viewing angle 
is measured with respect to the phase velocity, i.e. from a generatrix of the cone, and is {9 — a) if 
we measure 9 with respect to the physical velocity. As already pointed out in Paper I, arbitrarily 
fast variations can be observed when cosa = P and 9 = a. 

Oblique shocks with a ~ l/P are indeed expected in relativistic jets. If the external pressure 
drops suddenly, the jet decollimates and recollimates through a series of alternating conical shocks 
with this kind of aperture (Bowman 1994; Gomez et al. 1997). Furthermore, radio jets are observed 
to bend away from their initial direction already at the parsec scale; these deviations are probably 
amplified by projection effects, and are due to grazing collisions between the jet and the irregularities 
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of the ambient medium. In this context, one should note that the transverse pressure exerted by the 
jet on the ambient medium is made up of two parts, the comoving pressure and the ram pressure, 
of which only the latter increases strongly with F for incidence angles larger than l/F. It is then 
plausible that the external medium will give way until angles of this magnitude are reached. Of 
course, the shock will not be conical in this case; however, the variability is dominated by that part 
of the perturbation which is closest to the line of sight. All other parts add up on much longer 
timescales and contribute to the “steady” emission of the jet, so that a precise treatment of their 
timing is unnecessary. 

The condition that the viewing angle {9 — a) be small is not very restrictive; the sources which 
we study are selected a priori to belong to the blazar class, and the line of sight is selected a priori 
to be aligned with the jet at angles 9 <1/T (Urry & Padovani 1995); within the “blazing” cone, 
the major part of the solid angle resides close the cone surface, 9 ~ l/F, which is our requirement 
for fast variations. 

The length of the injection time determines the thickness of the slab representing the perturba¬ 
tion, and constitutes a firm lower limit for the variability timescale. It is tenable that the injection 
time is determined by the central engine, the black hole, and that the flow is modulated with the 
hole crossing time even at distances much larger than the hole radius. Here we follow an analogy 
with the models invoked for gamma-ray bursts (Rees & Meszaros 1992, 1994); the hyperrelativistic 
blast wave radiates the burst at a fraction of a light year, with a “global” timescale determined 
by the distance times 1/F^; on a finer scale, however, the light curve is modulated with a time 
typical of the central engine, a neutron star, whose dimensions have remained imprinted in the flow 
notwithstanding the enormous expansion. In the case of AGN’s, this minimum modulation time 
could be of a few minutes, compatible with the fastest variations observed. 

Fast variability at radio frequencies, in comparison with gamma-ray frequencies, has additional 
features which must be considered. The lifetime of the radio emitting electrons is much longer than 
the gamma-ray one; in the former case, then, at any given coordinate time radiation is emitted not 
only by the electrons which have just been shocked, but also by older electrons shocked previously, 
which (see Fig. 1) travel at different distances from the jet axis. Given the very long lifetime of 
the relevant electrons, the emission region can cover the entire cross section of the jet, even if one 
includes the decline of the emissivity due to the lateral expansion of the flow. The observed light 
curve is severely smoothed for a generic line of sight. In the following, we show that in a region of 
the parameter space fast variability is nonetheless possible, and that this region is ample enough 
to account for the observed statistics. 

A final caveat is related with the duty cycle of fast variability at radio frequencies. Intra Day 
Variables exhibit long periods of continuous activity, with many outbursts superposed at any given 
time, while at gamma-ray frequencies the outbursts seem isolated (observational constraints might 
play a role here). Our model is better suited to account for isolated outbursts, since the unfavored 
parts of each slab emit on longer timescales, and superpose with other slabs to raise the “steady” 
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level of the source. The fast varying emission is diluted to a small fraction of the total when many 
slabs are seen simultaneously. This is in qualitative agreement with the data, since the gamma-ray 
flux varies by large factors, whereas IDV is of the order of a few percent only. At any rate, a 
quantitative analysis requires a more realistic description than that of Paper I: here we include 
properly the dynamics of the flow across the shock, and the angular distribution of the radiation 
after the shock. 


3. THE FORMALISM OF THE MODEL 


The flow takes place within a funnel with rigid walls. The walls have an opening angle ijj with 
respect to the 2 axis, the symmetry axis of the system, and a relativistic beam of cold plasma 
impinges on them at an angle rj (Fig. ||). Before the shock the bulk velocity is ui ~ c and the 
Lorentz factor Ti >> 1. Since the bulk velocity after the shock, V 2 , must be parallel to the wall, 
the Rankine-Hugoniot conditions (see, e.g., Kennel &: Coroniti 1984 ) establish the Lorentz factor 
downstream, r 2 , and the position of the shock with respect to the initial direction of the beam - 
i.e. the angle a of Fig. |I]. In the fram e S* where the shock is perpendicular the Lorentz factor 
upstream is F^ = sin^ a + cos^ a, with F^ « Fi for angles a near 7r/2 and FJ < \/2 when 

a < 1/Fi. The shock is not relativistic in general, so the exact solution of the Rankine-Hugoniot 
conditions is complicated, and we approximate it as follows: 


* 

V2 = 

This expression reproduces the two limiting solutions: V 2 = v\/‘i ~ c/3 for a relativistic shock 
(TJ >> 1); and V 2 = v\/h in the opposite case (FJ < \/2), under the assumption that the electrons 
become relativistically hot even if the protons do not. In the laboratory frame, finally, the Lorentz 
factor downstream is: 

= I \ 2 — 

yl - (^)^cos2a-h (1 _ (a)2cos2a) 

where a is determined by making V 2 parallel to the wall: 


c I 41 1 16 r~ 

3y “ ^ (F*)2 ^ ^ (F*)4 


tan e = 



( 3 ) 


and e = a — r] is the angle between the shock and the wall. 

After the shock the electrons radiate effectively for a distance As < 2 ;. If t is the coordinate 
time, the plasma meets the shock in 2 ; = Vltcos{^jJ — e)/cos a, and the emission region is a ring 
in the x — y plane with internal radius r = vi t sin(^/ — e)/cos a. The non-zero particle lifetime 
gives the width of the ring : Ar = As (sin?/ — cos?/ tan(?/ — e)) (Fig. ||). The relation between 
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the coordinate time t and the observed time to ( the time at which a given wavefront reaches the 
observer) yields the infinitesimal surface element of the ring in t which contributes to the radiation 
observed in to'. 

(i'v 

AA{to,t) = Ay{to,t) \ — \t (4) 

dto 

Ay{to, t) = \l{r + Ar)2 — x = -^ [to — t(l — — —^cos0)] (5) 

* sm u c cos a 

The integral of Eq. ^ with respect to t ~ weighted with the radiation pattern P of Eq. ^ - gives 

the area Aw[to) of the source region that generates the flux in to, F{to)'-, this region is still ring- 

shaped, but lies in a plane inclined at an angle ^ with respect to z, tan^ = [(cos a/ cos{ijj — e) — 

^ cos0)/(^ sin0)]. The flux in to is proportional to A^ito), and we can describe the light curves 

at different frequencies by an appropriate choice of As. 


The high frequencies are radiated by energetic electrons, which cool very quickly so that 
Ar << r; then the model predicts sharp flares with timescales << At^ for a ~ l/Ei and for lines 
of sight close to a generatrix of the shock, 6 k. Op = if) — e (Paper I). At low frequencies the lifetime 
of the particles increases and Ar increases too, until it becomes comparable to r. The light curve 
is smoothed, and the shortest variability timescale approaches the crossing time Atz- In order to 
observe sharp temporal features the viewing angle 6 must be close to Op, and at the same time be 
such that ^{0) is close to V'- In this geometry As is parallel to the surface A{to) that contributes 
to the radiation at a given to', the photons emitted by the electrons during their lifetime reach the 
observer at the same instant, and the light curve is independent of As. The angle 0 for which this 
occurs is: 


cos a ± I tan ijjly (^)^ cos^ {tp — e)(tan^ tp + 1) — cos^ a 
^ ^^cos (V^-e)(l + tan2V:) 

Fig. § shows the theoretical light curves for different viewing angles around 0^. If we assume 
Ti = 15, a = l/Ti and 'ip = 0.1, then 0^ ~ 0.04, Op ~ 0.08, and r 2 ~ 10. If tq is the initial jet 
radius, the light curves are computed assuming that As = 10 tq: over this distance the jet radius 
doubles, and the emissivity drops by a factor of several. The variability is about 25% on a time 
scale ~ O.lAt^, if 0 is between 0.04 and 0.06. 


The light curves in Fig. ^ are computed with a realistic angular distribution of the radiation. In 
the frame comoving with the plasma the magnetic field upstream is assumed isotropic, downstream 
the magnetic held component parallel to the shock Bji dominates the normal one B±, even if 
TJ ~ 1. Including only By and transforming back to the laboratory frame, we get 


dP{x, 4>) 

(Kl 


1 + 


smecfcosy 


C ' 


1 - ^ cos X 


+ 


cos €c sm X 


r 2 (l - ^ cosx)\/l + tan^ (p^ 


r2 


( 7 ) 


where P is the power emitted in a unit solid angle in the direction dehned by x aiid (p, 6 = 
[r 2 (1 — V 2 /C cosx)] is the Doppler factor, Cc = tan(a — r])/T 2 is the angle between the shock and 














( 8 ) 


- 6 - 


the wall in the comoving frame after the shock, and 

X 

cos X = cos 9 cos V’(l H— tan 6 tan ijj) 


tan(j) = 


tan0-^l — (x/r)2 
sin ijj 


(9) 


If x',y',z' are cartesian coordinates with z' along V 2 (parallel to the wall), the z'x'^ and x'y' sections 
of the radiation pattern for ri=15 and a = l/hi are shown in Fig. If the opposite walls of the 
funnel diverge by an angle > the width of the radiation pattern, the long-lasting tail of a flare due 
to the ‘unfavored’ side of the funnel is much suppressed; this is important when modeling light 
curves which, as the radio ones, show repeated flares with a high duty cycle. 


4. COMPARISON WITH THE DATA AND DISCUSSION 


Fig. ^ shows the comparison of a theoretical light curve with the observations of S5 0716+714 


([Wagner et ah 19961) . This source is classified as a BL-Lac object, with a bright (S 5 ghz > 1 Jy) 
and flat radio emission [Su oc with a = 0.42 between 1.5 and 5 GHz, Ghisellini et ah 19971) , and 
a significant optical polarization. The optical spectrum is completely featureless and the redshift 
determination is difficult. The absence of any host galaxy rules out redshifts z < 0.25 and it is 
generally assumed z > 0.3. 

The light curve in Fig. Q represents the VLA observations of February 1990. The source was 
observed at A = 6 cm every two hours, with a regular sampling (except for weekly maintenance 
gaps). Normalized to the average brightness, the light curve shows intraday variability, in particular 
in the first part of the campaign. 

The theoretical light curve in the lower panel was obtained by convolving many curves like 
those of Fig. with random initial times and amplitudes. The time separation between successive 
microflares has a flat distribution between 0.025 and 0.075 xAt^, the amplitude has a flat top 
distribution of relative width 20%. The parameters of the elementary curves are equal to the those 
of Fig. with 9 = 0.04. The choice of the intrinsic time ro/(cFi) sets the physical size of the 
source and the variability timescale; here we have taken ro/(cFi)=13 days. Since similar results 
are obtained with larger 0’s up to 0.06, the visibility fraction is (0.06^ — 0.04^ ) X F 2 ~0.24, in 
good agreement with the occurrence of Intra Day Variability in the Effelsberg catalog. Indeed, our 
results are valid when the line of sight is at a small angle with respect to the shock surface; since 
the line of sight must be a priori within I/F 2 of the same surface, in order that the source shows 
blazar features, the visibility fraction is substantial. The model by Qian et al. (1991) is similar 
to ours in requiring an injection timescale shorter than the variability timescale; on top of this, 
the line of sight must be at a small angle with respect to the symmetry axis, instead of the shock 
surface, and the visibility fraction is diminished with respect to ours by a factor of order I/F 2 . 
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A final parameter to be chosen is the average flux, to which the curves are normalized. If the 
comoving brightness temperature is Tg, = 1/2 and Hq = 70 km s“^ Mpc“^, the flux emerging 
from the base of the cone in Fig. can be written as 

_ 2z/^(l + z)^A:Br2r^37r[cri 13 days/(l + 

(9 Gpc X (1 + 2: - Vl + z)Y 

The maximum value for Tg is 10^^ K, of course. At relatively low redshifts, appropriate to objects 
like S5 0716+714, we obtain fluxes of the order of those observed in the S5 catalog, i.e. about a 
Jansky. 

The particular source in question, S5 0716+714, was chosen as a paradigm of IDV at low 
redshifts because of the extensive data base, and from this point of view the performance of the 
model is comforting. However, S5 0716+714 is also the only firm example of correlated optical 
and radio IDV (WW). Most models have the high frequency emission coming from the inner jet, 
at distances of about 10^®-10^^ cm, and the radio emission from further out, at distances of about 
10^®-10^*^ cm. If the correlation is accounted for by making cospatial the optical and the radio, then 
the canonical scheme must be abandoned: placing the radio source in the inner jet would imply 
extreme brightness temperatures and coherent processes (see below). Moving the high frequency 
source further out would be compatible with the scheme proposed here, but would be difficult to 
reconcile with models involving inverse Compton, because of the low density of targets; one could 
perhaps resort to models involving proton induced cascades (Mannheim 1993). 

At redshifts of order 1 and higher, appropriate to quasar IDV’s, the fluxes predicted by Eq. 0 
are too low. We can obtain reasonable fluxes by taking higher values for ro/(cri), at the expense 
of longer variability timescales. These limitations amount to an upper bound on the apparent 
brightness temperature of several times 10^^ K (with the same definitions of WW). We conclude 
that, while a non-negligible fraction of IDV’s can in principle be accounted for by our model, the 
extreme events, and especially the 10^^ K flare observed in PKS0405-385, may involved alternative 
scenarios as coherence. 


This work was partly supported by the Italian Space Agency (ASI) throughout grants ARS- 
96-66 and ARS-98-116/22. 
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Fig. 1.— The source geometry. A beam of cold electrons (Fi) impinges on the jet walls at an angle 
r] . A shock is formed at an angle a. with the initial velocity of the particles, and the electrons 
became relativistic in the comoving frame (dashed part). They move along the jet walls (F 2 ) and 
radiate for a distance As. The viewing angle is 0. The ring in the figure is the emission region 
at the coordinate time t. Only the filled part of this ring (AA = Ay\dx/dto\t) contributes to the 
flux F{to)- The two dashed lines indicate the positions of the planes (normal to the page, and 
inclined by ^ with respect to z), from which the photons emitted at different t’s arrive together to 
the observer between to and to + dto- 
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Fig. 2.— The theoretical light curves for different viewing angles around 0^ = 0.042, with Fi = 15, 
a = 1/Fi, '0 = 0.1 and As = 10 tq. The time is normalized to the crossing time Atz and the ffux 
F{to) to its average value. If 6 is between 0.04 and 0.06, the variability is about 25% on a timescale 
of ss 0.lAt;J. 





Fig. 3.— The radiation pattern for Fi = 15 and a = 1/Fi. If x',y',z' are cartesian coordinates 
with z' along V 2 , the upper panel shows the section x'z'^ and the bottom one the section x'y'. 
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Fig. 4.— The observed (top panel) and theoretical (bottom panel) light curve of S5 0716+714, 
at 5 GHZ. The data are taken from Wagner et al. (1996). The abscissa is the time in days, and 
the ordinate the flux, normalized to the average one. The parameters of the theoretical curve are: 
Ti = 15, a = l/Ti, As = 10 ro, V' = O-l; ^’o/cFi =13 days. 





